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Abstract 



The equations of the mode-couphng theory (MCT) for ideal hquid-glass 
transitions are used for a discussion of the evolution of the density-fluctuation 



'^ ■ spectra of glass-forming systems for frequencies within the dynamical window 



between the band of high-frequency motion and the band of low-frequency- 
structural-relaxation processes. It is shown that the strong interaction be- 
T^ ■ tween density fluctuations with microscopic wave length and the arrested glass 

^^ ■ structure causes an anomalous-oscillation peak, which exhibits the properties 

^j. • of the so-called boson peak. It produces an elastic modulus which governs 

the hybridization of density fluctuations of mesoscopic wave length with the 
boson-peak oscillations. This leads to the existence of high-frequency sound 

^ ■ with properties as found by X-ray-scattering spectroscopy of glasses and 

O 

O ■ glassy liquids. The results of the theory are demonstrated for a model of 

- ^ . the hard-sphere system. It is also derived that certain schematic MCT mod- 

j^ ■ els, whose spectra for the stiff-glass states can be expressed by elementary 

formulas, provide reasonable approximations for the solutions of the general 

MCT equations. 

PACS numbers: 64.70.Pf, 63.50.+X, 61.20.Lc 



I. INTRODUCTION 

The excitation spectra of normal liquids are located in the THz band. But if the liquid 
can be supercooled new spectra evolve within the GHz band as precursors of the glass 
transition. These spectra are due to structural-relaxation processes. In addition there 
appear vibrational spectra for frequencies u above 0.1 THz which are characteristic for 
glasses and glass-forming liquids, namely high-frequency sound and the so-called boson 
peak. In this paper a microscopic theory shall be presented for these vibrational excitations. 

High-frequency sound was discovered only recently by X-ray-scattering spectroscopy in 
various glass-forming systems: in aqueous solutions of LiCl fTJI, glycerol [0,3, silica P-§, 
orthoterphenyl [|,g, 0.4 Ca(N03)2 0.6 KNO3 [|], and B2O3 [||. Molecular-dynamics simu- 
lations identified these modes for models of argon glass P,|lO| , ZnCl2 [|lll , and silica ||T2| , p!3 



Five features of high-frequency sound can be inferred from the cited work. First, the res- 
onance position Vt'f^^ of the dynamical structure factor Sq{uj) is a linear function of the 
wave number q^Vff^^ = VodQ, and this for q as large as 4 to 8 nm~^. Here v^o is the same 
high-frequency sound speed as known from Brillouin-scattering spectroscopy done for wave 
vectors which are about two orders of magnitude smaller. Second, the resonance width Tq 
exhibits a quadratic wave- vector dependence: Tq = 7^^. Hence density fluctuations propa- 
gate as known from the theory of elastic media for wave lengths down to the order of the 
inter-particle distances |jll|. This holds with the reservation that the spectrometer resolution 
of about 0.4 THz causes considerable uncertainties in the data for f2™^^ and Tq and that it 
is difficult to separate the inelastic-scattering signal from the huge background due to the 
resolution-broadened quasi-elastic scattering. The damping of ordinary sound in crystals or 
liquids depends strongly on temperature. As a third remarkable feature of high-frequency 
sound one observes that Tq exhibits only a weak temperature dependence. For silica and 
orthoterphenyl no temperature dependence of Tq could be detected even though fi™'^'' de- 
creases with increasing temperature considerably 0,^,|^,|^. For hydrodynamic sound the 
resonances exhaust the inelastic spectrum of the density fluctuations. Contrary to this one 



infers from the simulation work [^,^,0 as a fourth property of high-frequency sound that 
the resonances are superimposed on a broad intense background spectrum of Sg{uj). As a 
fifth property one finds for the simulation results of silica |13[ that the background spectrum 
exhibits a low-frequency threshold. This conclusion is supported by the measurement of the 
X-ray-scattering spectrum of densified silica j^, which demonstrates that there is a sup- 
pression of Sq{uj) for small u relative to the hydrodynamics description for the dynamical 
structure factor. 

The so-called boson peak has been known for quite some time since it can be detected 
in the spectra obtained by standard techniques. It was measured for all the above cited 



systems, for example by Raman scattering for LiCl [|I4| and by neutron scattering for silica 
fljl5[| . It was found by molecular-dynamics simulations for ZnCl2 |1^ and silica |lT7| , p!8| . Six 
features are typical for these peaks. First of all, the peaks are due to soft excitations. The 
position Up of the peak maximum is several times smaller than the Debye frequency ujd 
of the system. In silica the peak causes an enhancement of the so-called density of states 
by a factor of about seven relative to the Debye spectrum |TH]. Second, the anomalous 
peak is due to quasi-harmonic oscillations. Originally, this was concluded indirectly, since 
this assumption can explain the enhancement of the specific heat above Debye's T^-law 



for temperatures T near lOK [1^. Normal-mode analysis for a model of silica showed 
this result explicitly [0. Third, the peak is skewed. There seems to be a low-frequency 
threshold i7_ so that the high-frequency wing of the peak extends further than the one on 
the low-frequency side [0,0,|1^]. Fourth, the peak appears to be related to some instability 
of the system since up and r2_ decrease to zero upon heating the system towards some 



characteristic temperature Ti,p [|1^,0. The mentioned four features are also exhibited by 
the density of states of a harmonic lattice with a random distribution of force constants 
p!9| . As a fifth property one observes that a quasi-elastic relaxation peak appears for T 
increasing towards Tbp which eventually buries the peak |T^ . Sixth, the peak position ujp of 
the dynamical structure factor depends only weakly on wave vector q if at all [§|,|TB . 



The discussions of this paper are done within the framework of the mode-coupling theory 



(MCT) for the evolution of glassy dynamics in simple systems. This theory is based on 
closed microscopic equations of motion for the density-correlation functions pO| , |2T| . Previous 
work done in this context focused on elucidating the properties of the structural-relaxation 
phenomena, as described in Refs. p^j23|| and the papers quoted there. In the following it 
will be shown that MCT also implies a theory for the evolution of anomalous-oscillation 
peaks (AOP) with the same six properties specified in the preceding paragraph for the so- 
called boson peak. Furthermore it will be demonstrated that the AOP manifests itself in 
resonances of the density-fluctuation spectra for wave-vectors up to half the value of the 
Debye vector qu-, which exhibit the five features listed above for high-frequency sound. 

MCT deals with states of matter for which the structure factor Sq depends smoothly 
on wave vector q and on control parameters like the temperature T and density p. The 
equations of motion of the basic version of this theory, which will be used in this paper, 
exhibit a bifurcation from ergodic liquid dynamics to non-ergodic glass dynamics if T or 
p cross critical values T^ or p^ respectively. This bifurcation provides a model for an ideal 
liquid-glass transition. In the extended version of MCT the singular transition is replaced by 
regular crossover PB[ . The crossover is connected with the evolution of structural relaxation. 



Extensive tests of the MCT predictions have been performed in recent years. Let us only 
mention one group of such tests which are of particular relevance for the following discussions. 
The MCT results for the Debye- Waller factors fq of the glass depend non trivially on q and 
on T or p. The results calculated for the hard-sphere system [^ agree with the data 



measured for hard-sphere colloids [24]. Similarly, the /„ calculated for a Lennard- Jones 



mixture p5[ agree within 10% with the values deduced from molecular-dynamics studies 
26| . These assessments of the theory and many other tests, which are reviewed in Ref. 
27| , show that MCT describes properly some essential features of structural relaxation. 
Since the MCT equations are rather involved most of the original work was done first for 
schematic models. These are truncations of the large set of equations to sets dealing with a 
few correlators only. The truncations are done with the intention to reduce the complexity 
of the mathematical machinery without loosing the essence of the features to be analyzed 



for the solution. It was noticed for schematic models dealing with a single correlator that 
the spectra for the states for T <^ T^. or p ^ pc exhibit a broad peak, which is due to a 



superposition of harmonic oscillator spectra [22,E8l. Tao et al. recognized that the evolution 



of these peaks due to changes of control parameters was similar to what they had measured 



for the boson peak of glassy aqueous LiCl |T^. Therefore they concluded that MCT implies 



a theory for the boson-peak spectrum. This conclusion was corroborated by a series of 
studies, where solutions of schematic models were used to describe quantitatively spectra 



of glassy liquids p9H3^. The fits described structural-relaxation spectra in windows of 
several orders of magnitude in size in addition to parts of the boson-peak spectra for high 
frequencies. In the following this earlier work will be extended to a systematic theory for 
the general MCT equations. This theory implies, in particular, the suggestion that the 
characteristic temperature Tf,p for the boson-peak dynamics is identical with Tf.. Our results 
will be demonstrated comprehensively for a model of the hard-sphere system (HSS). 

The paper is organized as follows. In Sec. II the basic equations for our calculations 
are formulated and the details for the HSS work are specified. The evolution of structural 
relaxation will be demonstrated in order to put the discussion in the proper context with 
the theory of the glass transition. Then, in Sec. Ill, the results for the evolution of the 
anomalous-oscillation peak and for high-frequency sound are presented. These oscillation 
features of the MCT dynamics are described within a generalized-hydrodynamics approach. 
In Sec. IV the phenomena are explained within a schematic-model analysis. In the concluding 
Sec.V our findings are discussed. 

II. BASIC EQUATIONS 

A. A mode coupling theory model for a dense simple system 

The basic quantity specifying the equihbrium structure of a simple system is the structure 
factor Sq = {\ pg- \'^), which is the canonical average of the squared modulus of the density 



fluctuations of wave vector q : p^ = Y.a^w{iQ^a)/VN. Here a labels the N particles at 
positions fa in the system of density p. The structure factor of amorphous systems depends 
on the wave vector modulus q =\ q\ only and it is usually expressed in terms of the direct 
correlation function Cg : Sg = 1/(1 — pCg) ||3^. The most relevant variables describing the 
dynamics of structure changes as function of time t are the density correlators (f)q{t) = 
{p*^(t)Pq) / Sq. The short time asymptote of these functions is specified by a characteristic 
frequency Vtq : (pqit) = 1 — i(i7gt)^ + ■ ■ ■. In the small-wave-vector limit one gets the 
dispersion law for sound; Vtq = Voq + O(g^), where Vq denotes the isothermal sound speed. 
For general wave vectors one obtains Q^ = v'^q'^/Sq, where v denotes the thermal velocity of 
the particles |3^. In particular Vq = v'^/Sq. Within the Zwanzig-Mori formalism one can 
derive the exact equation of motion 

0,(t) + fij0,(t) + nl f ruqit - t')<Pq{t')dt' = . (la) 



The relaxation kernel rriqii) is a correlation function of fluctuating forces |^. Let us intro- 
duce here and in the following Fourier-Laplace transformations to map (f)q,mq and similar 
functions from the time domain onto the frequency domain according to the convention 
(pq{uj) = i J^ exp{iu!t)(f)g(t)dt = (p'qioj) + i4>g{uj). The reactive part 4>'g{uj) and the dissipa- 
tive part or the spectrum (f)'g{uj) are connected via a Kramers-Kronig relation. Equation 
(p!a|) is equivalent to the representation (pqi^uj) = —l/[uj — ^l/[uJ + Qgrrig^uj)]]. The fluctua- 
tion dissipation theorem connects (t>q{uj) with the dynamical susceptibility Xgi^) '■ Xgi^) = 
[u!(f)q{uj) + l]x^- Here xj = Sq/{ppv'^), with p denoting the mass of the particles, is the 
isothermal compressibility [^]. Therefore Eq. (|T^ is equivalent to 

xM/x; = -^v [^' - fi? + fi;^m,(c.)] . (lb) 

Within the MCT kernel mq{t) is written as the sum of a regular contribution and a 
mode-coupling contribution, describing the cage effect of dense systems within Kawasaki's 
factorization approximation. In this paper we neglect the regular term. Hence one gets ||2U[| : 



m,(t)=^,[0(t)] , (2a) 

where the mode-couphng functional jFg, considered as functional of the dummy variable /, 
is: J^q[f] = J2k+^=q'^iQy^P)fkfp- The coupling coefficients V{q,k'p) are determined by the 
equilibrium structure: 

V{q, kp) = pSgSkSp [q{kck + pCp)J /(2g ) . (2b) 

Let us approximate wave-vector integrals by Riemann sums, obtained by choosing an 
equally-spaced wave-vector grid of M terms. Thus the wave-vector index q will be un- 
derstood as a label running from 1 to M. Correlators 0(t), kernels etc. are considered as 
M-component vectors. After this discretization the functional is a quadratic polynomial 

M 

^,[f] = E Vg,kjkfp . (2c) 

kp=l 

The positive coefficients Vg^kp are related trivially to V^qjkp) in Eq. (^ ); for details the 
reader is referred to Ref. p6| . 

Equations (|l]) and (|^) are closed. They define a unique solution for all parameters 
Qq > 0, Vq^kp > 0. The solution has all the properties of a correlator, i.e. the (f)q(t) are 
positive-definite functions given by non-negative spectra (pg{uj). And the solutions depend 
smoothly on Qq and V^,fcp for all finite time intervals |37,38|. Thus, Eqs. (|l]) and (|^) formulate 



a well-defined model for a dynamics. The model is related to the physics of liquids by 
specifying Qg and Vq^kp in terms of particle interactions and the control parameters p, T. 

In the following the results will be demonstrated for the hard-sphere system (HSS). 
The structure factor is independent of T in this case, so that the only non-trivial control 
parameter is the packing fraction (f of the spheres of diameter d : (f = nd^p/G. The 
structure factor Sq is evaluated in the Percus-Yevick approximation [^. Wave vectors will 



be considered up to a cutoff value q* = AO/d and they are discretized to M = 300 values. 
The sphere diameter will be chosen as the unit of length, d = 1, and the unit of time will 
be chosen such that the thermal velocity is w = 2.5. 



B. Ideal glass states 



The specified model exhibits a fold bifurcation [^. For small coupling constants V^^p 



the correlators (f)q{t) and kernels mq{t) decay sufficiently fast to zero for times tending to 
infinity so that the spectra 0''(^) and m''{uj) are continuous in uj. Density fluctuations, 
which are created at time t = 0, disappear for long times; and the same holds for the force 
fluctuations. The system approaches equilibrium for long times as expected for an ergodic 
liquid. This implies lim^^Q ujiriqluj) = and one concludes from Eq. (^) that the static 
susceptibility Xq = Xq{^ -^ 0) agrees with the thermodynamic susceptibility, Xq = x'^ ■ For 
large coupling constants, on the other hand, there is arrest of density fluctuations for long 
times: (pqit -^ oo) = fg] < fg < 1. Thus the perturbed system does not return to the 
equilibrium state. Similarly, there is arrest of the force fluctuations: mg{t -^ oo) = Cg > 0. 



The numbers fq,Cq are connected by the equations pOf : 

fg = Cq/{l+Cq) , Cq=J'q[f] , g=l,2,...,M. (3) 

For this strong coupling solution the kernel exhibits a zero-frequency pole, limi^_>o ujmq{uj) = 
—Cq, and therefore one concludes from Eq. ([TED that the static susceptibility is smaller 
than the thermodynamic one, Xg < xj; since Xq = xj/(l + Cq). This is a signature for 
a non-ergodic state |39|. The system reacts more stiffly than expected for a canonical 
averaging. The dynamical structure factor Sq{uj) = Sq(t>'g{uj) exhibits a strictly elastic peak: 
Sq{uj) = -KSqfq5{uj)+ rcgular terms. This is the signature for a solid with fq denoting its 
Debye-Waller factor. Hence, the strong coupling solution deals with a disordered solid; it is 
a model for an ideal glass state. If one increases the coupling constants smoothly from small 
to large values one finds a singular change of the solution from the ergodic liquid to the 
non-ergodic glass state, i.e. an idealized liquid-glass transition. For simple-liquid models 
the transition occurs upon cooling at some critical temperature T^ or upon compression at 
some critical packing fraction ip^ pO|. For the HSS model under study one finds (pc ~ 0.516 
]. \i Lf increases above ip'^, the Debye-Waller factor fg increases above its value at the 



critical point, called the plateau /^, as is demonstrated in Fig. |I[ 
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In the theory of crystalhne sohds one defines susceptibihties with respect to the restricted 
ensemble of a given arrested lattice. The MCT equations of motion allow for a similar 



formulation of the equations of motion for the glass state |Q. To see this, one has to map 
the density correlators (f)q to new ones (j)q by 

Ut) = /, + (!- U)Ut) ■ (4a) 

If one introduces new characteristic frequencies Clg by 

^l = fi?/(l - /.) , (4b) 

one obtains the short time expansion (pq{t) = 1 — |(r2gt)^ + ■ ■ ■ in analogy to what was found 
for 0g(t). Substitution of these results into Eq. (0) reproduces the MCT equations of 
motion with (pg, Qg and m^ replaced by 0g, fig and m^, respectively. Here the new relaxation 
kernels are related to the original ones by 

mg{t) = C, + (1 + Cg)mg{t) . (4c) 

For the dynamical susceptibility one obtains the formula Xq{^) = Xq[l +^4>q{^)]- The new 
correlator has a vanishing long time limit, 0g(t —>■ oo) = 0, the Fourier-Laplace transform 
exhibits the property lim^^QUcpg^uj) = 0. Since the equation of motion does not change its 
form, one gets as the analogue of Eq. ([TEI ) the expression of the susceptibility Xq{^) iii terms 
of the polarization operator rhq{uj): 

Xq{^)IXq = -^V [^' - ^; + ^l^mg{uj)\ , Xg = x^l - fq) ■ (4d) 

Combining Eq. (^ with Eq. (^ one concludes rhg{t -^ oo) = 0, i.e. the new kernel rhg{uj) 
for the amorphous solid exhibits a regular zero-frequency behaviour: lim^^o^^q('^) = 0. 

The mentioned equation of motion for 0g(t) can be closed by combining Eqs. (^a]) and 
(^ with Eqs. (^) and (|2^). One finds an expression of the new kernel as a new mode- 
coupling functional JF of the new correlators: 

rhg{t) = fg U(t)l . (5a) 



j^ is a. sum of a linear term JF(^) and a quadratic one JF^^), JF^ = j^j-^^ + ^j;'^\ where: 

^!>'\f] = J2v,,kfk , v,,k = 2ii-f,)j2v,,kpm-fk), (5b) 

k P 

•^f [/] = E ^.,afc/p , ^.,/^P = (1 - f,)y<i,kpii - fk)ii - u) . (5c) 

kp 

As a result equations of motion are produced, which are of the same form as the MCT 
equations discussed in Sec.IIA. But in addition to the quadratic mode-coupling term there 
appears a linear one, and the values of the mode-coupling coefficients from Eq. ( pE| ) are 
renormalized to the coefficients Vq^kp- 

The preceding Eqs. (|]) are equivalent to equations relating the density-fluctuation 
spectra 0''(co') with the kernel spectra rh''{uj) = rh^^^"{uj) + m^^^"(uj): 



mi}^"(to) 



1 

k 



EWfc(^), (5d) 



kp 



-J2 duji duj2Vq,kpS{uj -uji- uj2)(j)l{uJi)(i)p{uj2) ■ (5e) 



One can interprete Eqs. @) and (H) using the language of the theory of boson fields, e.g. 
the phonon fields in crystals. Xq{ijj)/xq is the field propagator and Vtq is the bare-phonon- 
dispersion law. The kernel rhq{uj) is the phonon self energy. Equations ( pd|) and (|5^ ) are 
golden-rule expressions for the phonon-decay rates. Kernel m^^^ describes elastic scattering 
of the phonon from the disorder, produced by the amorphous glass structure. Kernel m*^^-* 
deals with the decay of a phonon into two due to anharmonicities. The glass structure 
infiuences the decay rates via the Debye-Waller factors which enter Eqs. (^) and (^ 
for Vq^k and Vq^kp, respectively. The challenge is to evaluate this probability rh'q{uj) self 
consistently; the decay depends on the same phonons as one wants to study. 

C. The glass— transition scenario 

Figure ^ exhibits the evolution of the dynamics with increase of the packing fraction 
(f. The wave vector q = 7.0, used in the lower panel, is close to the structure-factor-peak 
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position where the static susceptibihty x^ ^ Sg is high, while the wave vector q = 3.4, used 
in the upper panel, deals with fluctuations where x^ is very small (compare Fig. [|). Figure 
^ shows the equivalent information for the fluctuation spectra <f)''{uj). In Ref. |^ a further 
set of diagrams for the wave vector q = 10.6 can be found. 

The curves for ip < ipc with label n = 1 in Figs. ^ and ^ refer to the packing fraction 
(f = 0.276. The correlators exhibit strongly-damped oscillations; and the ideal phonon 
resonances, to be expected for nig = 0, are altered to broadened bumps in the spectra. If 
the packing fraction increases into the interval 0.9 fc'£f< fc, the oscillatory features 
almost disappear for t > 0.2; and the shown spectra (pgiuj) decrease monotonously with 
increasing uj. Simultaneously, the decay to equilibrium is delayed to larger times. At the 
critical point the correlators approach the plateau f^ in a stretched manner, as shown by the 
curves with label c in Fig. 0. This process, which is called critical decay, leads to a strong 
increase of the fluctuation spectra with decreasing frequency, as shown in Fig. |^. Increasing 
if above the critical packing fraction ip^ the values for the long time limits fg increase. Since 
these limits are approached exponentially fast for ip ^ ip^ [^], the correlation spectra become 
cij-independent for low frequencies. The value for this white noise spectrum increases if cp 
decreases towards (pc] this is a precursor phenomenon of the glass melting at the transition 
point. If (p is sufficiently far above (pc, one observes oscillations again, as can be seen for 
the ip = 0.6 result in the upper panel of Fig. |^. In this case, the correlation spectra are not 
anymore monotone functions of frequency. For q = 3.4, there occur two peaks for uj > 10. 
The narrow peak is due to high-frequency phonon propagation. In addition there is an 
anomalous-oscillation peak (AOP) for u ~ 80. For g = 7 a phonon peak is absent in the 
spectrum for (p = 0.6 but an AOP is present, as is shown in the lower panel of Fig. ^. 

The time scale for normal-state-liquid dynamics is set by the Debye frequency un- It 
is the same scale as for the dynamics of the crystalline state of matter. For the discussion 
of this normal condensed-matter dynamics it is sufficient to consider a window of, say, two 
decades around 271 /uc- This is demonstrated in Fig. g for the liquid state with label n = 1 
or for the glass state with label n = 3. If the time increases from 0.01 to 1 the correlators 
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(f)q{t) decrease from 0.9 to the long time limit. All oscillations occur within this interval, 
which is also called the transient regime. The spectra for these states are located within 
a corresponding regime of microscopic excitations, extending roughly between 0.3 and 300. 
For frequencies around and below 0.3 there is only white noise for the specified normal 
states. The glass transition is connected with a dynamics, called glassy dynamics, which 
occurs for times longer and frequencies smaller than the ones characterizing the transient. 



For reasons which can be understood by asymptotic expansions about the critical point p2 
the dynamics is stretched over many decades. The glassy dynamics of the HSS is discussed 
comprehensively in Refs. [^,^. It is impossible to view glassy dynamics adequately on 
linear scales for t or u. Conventionally, one represents the results on logarithmic abscissas 
as done in Figs. ^ and ^ The bifurcation for (f = (f^ or T = T^ also modifies the transient 
dynamics. The subject of the paper is the study of these modifications. 

III. EVOLUTION OF THE TRANSIENT DYNAMICS 

A. Anomalous oscillation peaks and high frequency phonons 

For a discussion of the transient motion it is sufficient to consider dynamical windows 
of about two orders of magnitude. These can be viewed more adequately on linear rather 
than on logarithmic abscissas. There is no reason to consider such fine tuning of control 
parameters relative to the critical point as is necessary for a study of structural relaxation. 
Therefore let us extract the relevant information for the following discussion from Figs. |^ 
and |] and replot it as Figs. § and ^. Let us first consider the results for wave vector q = 7. 
For the packing fraction ip = 0.60 the particles are localized in such tight cages that the 



square root Sr = J {5r^{t -^ oo)) of the long-time limit of the mean-squared displacement 
is only 5.0% of the particle diameter [^. Thus one expects the particles to bounce in 
their cages with an average frequency of the order u ^ 2'kv/{A ■ 6r) = 78. This explains 
qualitatively the oscillation around the equilibrium value /7.0 exhibited by 4>7.o(t) in Fig. ^ 
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for t < 0.3 and the corresponding "peak" of the spectrum in Fig. |^. The maximum position 
ujp K, 75 of this AOP is estimated well by the crude formula. The AOP differs qualitatively 
from the Lorentzian which one would expect for some damped harmonic oscillation. The low 
frequency part of the peak decreases more steeply with decreasing uj so that there appears 
some threshold near uo = 40. Below the threshold the spectrum is flatter than expected for 
the wing of a Lorentzian. If (f decreases to 0.567 the cages widen so that 6r = 0.070 |^ . 



This explains the increase of the oscillation period exhibited by curve n = 3 in Fig. ^ and the 
corresponding downward shift of up in Fig. ^. The shift is accompanied by a strong increase 
of the spectrum for uj < 10. The integral of the inelastic spectrum also increases, reflecting 
the decrease of the elastic contribution vr/^ which is exhibited in Fig. |l|. The described 
trends continue if (p is decreased further to 0.540 (curves n = 4). The threshold and the 
spectral minimum for small frequencies are now replaced by a central peak for cj < 20. The 
maximum of the AOP, estimated from 6r = 0.099 [^ as up ~ 40, is buried under the tail 



of the central peak; it merely shows up as a shoulder. The correlator still exhibits some 
small oscillation before it reaches its long time limit /7.0 = 0.96 for t > 0.2, but it does not 
fall below /7.0 anymore. At the critical point 6r = 6r'^ = 0.183 |^ and so one estimates a 



position Up ?^ 21. But the critical decay manifests itself by the appearance of a long time 
tail of 07.o(^)- The approach to the asymptote /70 = 0.85 cannot be demonstrated on the 
linear time axis used in Fig. ^. This tail leads to a strong enhancement of the central peak, 
so that the AOP cannot be identified. This trend continues if a packing fraction below the 
critical value is considered. 

For the glass states the correlators 03.4 (t) exhibit weakly damped oscillations which 
lead to nearly Lorentzian resonances for the spectra. These excitations are analogous to 
phonons in crystals. The softening of the glass with decreasing ip leads to a decrease of 
the phonon frequency, but Fig. ^ demonstrates that for all p > Pc the peak positions 
are located considerably above the maximum position ujp of the AOP discussed in the 
preceding paragraph. One recognizes for the p = 0.60 result in Fig. ^ that the oscillations 
for t < 0.1 do not occur around the long time limit /a, 4. Rather the oscillation centre follows 
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a curve discussed above for the bouncing in the cage. This is equivalent to the fact that 
the phonon resonance does not exhaust the spectrum (f)'^^{u), rather it is placed on top of 
some background. The background exhibits a similar threshold as discussed above for the 
AOP of 07 (a;), and with decreasing ip it also follows the same pattern as described for the 
q = 7.0 spectra. Apparently, the dynamics for q = 3.4 illustrates a hybridization of the 
phonon dynamics with the modes building the AOP. The regular dependence of the MCT 
solutions on control-parameter variations are the reason why the results in Figs. § and 
P do not exhibit any drastic change if (f is shifted from the glass state with label n = 4, 
referring to e = [if — ipc)/'^c = 0.0464, through the critical point to the liquid state with 
e = —0.0464. But upon shifting the state into the liquid the phonon resonances get buried 
under the relaxation spectra. 

The variation of the spectra with changes of the wave vector q is demonstrated in Fig. ^ 
for the stiff-glass state with the packing fraction (f = 0.60. For q < 0.6 a single peak of nearly 
Lorentzian shape exhausts the whole inelastic spectrum (pgiuj). The resonance position 
follows the dispersion law of high-frequency sound Q^'^'^ = Vooq. Here t^oo is the sound speed 



expected from the glass susceptibility Xg=05 Eq- P^, ^oo = Vo/Jl — fq=o = 75.8. Also the 
half width of the resonance exhibits the quadratic wave-vector variation expected for sound 
in an elastic continuum: Fg = 7 ■ g^. The single-peak shape of the spectra is found for 
all wave vectors up to about qof^, where q^ = {367r(fY^'^ = 4.08 denotes the Debye wave 
vector of the system. Also the linear dispersion law continues up to these large g-values as 
is demonstrated in Fig. |^. However, for q > 0.6 the resonance width Tg is somewhat larger 
than expected by extrapolating the asymptotic law F^ = '^{Q^^^/vooY, as is shown in Fig. 
p. The sound frequency Q"^^^ reaches the position of the maximum of AOP at g ?^ 1.2 and 
extends considerably beyond this value for larger q, as is shown for g = 1.8 in Fig. |^. 

If fi™^^ increases to the centre of the AOP the loffe-Regel limit is approached; i.e. the 
sound frequency becomes of the same order as the resonance width. In this case hybridization 
of high-frequency sound with the modes forming the AOP becomes important. The sound 
resonance is not Lorentzian anymore, as demonstrated in Fig. ^. For q = 1.4 the threshold 
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of the AOP near u = 40 causes a shoulder on the low-frequency wing of the sound peak. 
With increasing q a broad flat background spectrum is formed between the threshold and 
Q^^^ as shown in Fig. |^ for q = 1.8 and 2.2. For wave vectors exceeding g_D/2 the spectrum 
exhibits two maxima. The one at lower frequency is due to the AOP. The peak at higher 
frequencies is the continuation of the high-frequency sound resonance, which is increased 
because of a level repulsion effect, as demonstrated in Fig. ^ for q between 2.2 and 5.4. 
The high-frequency-phonon frequency Vt^'"^^ increases with q increasing up to about qfj, 
and then it decreases again if q increases up to about l.^qn. So the dispersion follows the 
pattern expected for a crystal phonon near the Brillouin zone boundary. The width Tg of 
the high-frequency phonon for 2.5 < q < 5.0 varies only weakly with changes of q. Also the 
AOP maximum ujp does not change very much if q increases from 2 to 5, as is demonstrated 
by the open circles in Fig. ^ In a large wave-vector interval around the position of the first 
sharp diffraction peak of the structure factor Sg, the spectra do not exhibit a phonon peak 
but the AOP only. The spectral shape varies only weakly with changes of q, except near 
g = 8.7; there a new splitting in two peaks occurs. 

In Fig. ^ the spectroscopic parameters Q"^^^ and Tg of high-frequency sound are pre- 
sented as function of the packing fraction. The resonance frequency decreases with decreas- 
ing if, reflecting the softening of the glass state upon expansion. Remarkably, the damping 
Tg does not vary much with changes of cp. 

B. A generalized hydrodynamics description for the glass state dynamics 

The description of hydrodynamic sound is obtained by coarse graining correlation func- 
tions so that fluctuations on microscopic scales for the space-and-time variations are av- 
eraged out. Coarse graining is equivalent to restricting functions like Clg or rhg{uj) to their 
leading order Taylor coefficient with respect to their g-and-tu dependence. In this manner 
one gets from Eq. d^d]) the hydrodynamics description: 

xU^)lx, = fif V h' - fif + ^^rf ] . (6a) 
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Here the dispersion law Vt^ and the damping T^ are given by 



^f = ^^oo? 



vl{l + C) ,rf = 7g' , i = vlK'\u: = Q) 



(6b) 



where the notation C = Cq=o and K{t) = rhg^oit) is introduced. Because of Eqs. (^), (^ 
and (^) one can express C and -ft'(t) in terms of the mode-couphng functional jFg = J^q=o: 



C = ^o[/] , K{t) = { J-o[0(t)] - C} /(I + C) . 



An explicit expression for the new functional JFq follows from Eqs. 



•^o[/] = / V,f,dk 
Jo 



(7) 



and dbj): 



(8a) 



where the weight factors V^ > read 



Vk = pSg=o [Skk/2'K] 



cl + ^{K)ck + -{K)^ 



(8b) 



The integral in Eq. (|8^) is to be understood as a Riemann sum over the wave-vector grid 
of M terms as specified in Sec.IIA. 

The hydrodynamics description is not suited to deal with high-frequency sound and 
the AOP. In order to identify the essence of these phenomena a generalized-hydrodynamics 
description shall be developed. It is obtained from Eq. ( ^d]) via approximating kernel rhq{uj) 
by its zero- wave-vector limit: 



xr''M/x. = ^;/ co'-nlil-coKico)) 



(9) 



This formula extends the hydrodynamics approximation in two respects. Firstly and most 
importantly, the hydrodynamic damping constant i'y/v'^ is replaced by the frequency- 
dependent function K{uj) = K'{uj) + iK"{uj). Secondly, the full non-trivial q dependence 
of Clq is kept. According to Fig. |^ it is the g-dependence of Clg which dominates that of 
the whole density-fluctuation spectra. If one intends to describe these spectra for wave 
vectors up to the inverse of the interparticle distances, one must not simplify Clq. However, 
if one is merely interested in the description of high-frequency sound so that q is restricted 
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to the regime below qo/'^i say, one can replace Ctq in Eq. (P) by the hydrodynamics limit 
^q = v^q. In this case, formula (H) has a form often used in acoustics, where {1 — ujK{uj)) is 
called the dimensionless longitudinal elastic modulus. Let us also rewrite the mode-coupling 
formula (0) for K{t) more transparently in analogy to Eqs. (j^): 

K{t) = K^^\t) + K^^\t) , (10a) 

K(i)(t) = r W^«4(t) , 1^« = 2(1 - /,=o)^./fc(l - fk) , (10b) 



i^(2)(t) = r W^f 0,(t)2 , ^f = (1 - fq^,)V,{l - hf . (10c) 

^ u 

In Fig. 1^ representative MCT spectra for the glass states of the HSS are compared 
with the results of the generalized-hydrodynamics approximation. For the regime q <2 the 
description by Eq. (^ with Clg = v^oq is very good. This does not only hold for the treatment 
of the high-frequency-sound resonance but also for the background spectrum. But the 
generalized-hydrodynamics approximation yields reasonable results also for wave vectors 
extending up to and exceeding the value of the structure-factor-peak position g^ax ~ 7. 
In particular the subtle hybridization of the phonon with the modes underlying the AOP 



are treated semi quantitatively correctly. Hence, one concludes that Eqs. (P) and (|ToD 
separate the problem of the evolution of the transient dynamics into two parts. The first 
part concerns the dependence on the wave vectors. This is determined entirely by the 
quantities Xg and Qg. These functions are constructed from the structure factor Sg, whose 
dependence on q and ip is well understood. Sg determines the characteristic frequency fig and 
the thermodynamic susceptibility x^- The MCT equations provide the glass form factors 
fg, which are understood as well [^. And the combination of this information yields Clg 
and Xg via Eqs. ( ^ ) and ( ^ ) . Consequently, the remaining issue of this paper is to provide 
the understanding of the kernel K{uj). 

Figure 110 exhibits the spectra K"{uj) and also the reactive parts of the modulus A(u;) = 



1 — ujK'{uj) for three representative glass states. The spectra K"{uj) are much simpler 
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than those for the correlators since phonon resonances are absent. This property is the 
essential result achieved by the Zwanzig-Mori theory, and MCT preserves this property. 
The spectra of the kernels consist of some broad background and an AOP. The AOP of 
K"{uj) has quite a similar form as discussed in Sec.IIIA for (J)"q{uj) and also its changes due 
to changes of (p are similar. Formula (^ describes the hybridization of two "oscillators". 
One is a bare phonon with dispersion Clg and the other is represented by the AOP spectrum. 
The hybridization problem is analogous to the one considered, for example, in the dynamical 
theory of light propagation in dielectric media. The bare phonon corresponds to the electro- 
magnetic wave in vacuum and 1/[1 — ujK{uj)] is the analogue of the dielectric function. 
In lowest order one gets two resonances of the coupled system whose frequencies Q'^'^^ are 
obtained approximately as solution of the equation (fi^^^/fig)^ = A{Cl^^^) and the damping 
is estimated by Fg = Q'^K" (Q^^^^) . Elementary discussion of these equations with the aid of 



Fig. |T0|, which is left to the reader, explains qualitatively the findings reported in Figs. 

A digression might be useful concerning the different role played by the two equations 
for the susceptibility ([T^) and (|15D. The crucial property of glassy dynamics is that kernel 
mg{uj) in Eq.(|r^) depends strongly on frequency u. In particular, this kernel has a large 
reactive part m' Ju). Therefore it makes no sense to try a hydrodynamics approximation for 
Eq. ([Tbl) , based on nig^uj) ^ niq^uj = 0) = im'g^u = 0). But let us consider a generalized- 
hydrodynamics approximation, defined by ignoring the wave vector dependence of kernel 



mq{uo). Because of Eq. (|2a|) this limit is given by jFo[0(t)] and thus it can be expressed by 
K{uj) via Eq. ©: 

mg=o(^) = [-C/uj] + (1 + C)K{uj) . (11) 

The kernel mg=o{u) exhibits the pole [—C/u], which is the signature of the ideal glass state; 
other subtleties of the glassy dynamics are hidden in K{uj). Substitution of Eq. (|Tl]) for 
mq{uj) in Eq. ( [151) reproduces Eq. (|^), except that the renormalized frequency Clg is replaced 



by the approximation to Eq. (^): Vtq/ Jl — fq=o. For our model of the HSS structure there 
are no serious wave- vector dependencies of fg for q < qof^- This means that high-frequency 



sound and the background spectrum for g < 2 can also be discussed on the basis of Eq. (|T^). 
The dotted hues in Fig. ^ exhibit the corresponding resuhs for q = 1.0 and 3.4. However, 
Fig. |I] demonstrates that the relevant factor l — fq = 1/(1 + Cq) varies by more than 100% for 
q approaching and exceeding qo- Therefore the generalized-hydrodynamics approximation 
based on Eq. ( ]TT| ) cannot be used to discuss, for example, (pjQ^uj). The superiority of the 
generalized-hydrodynamics approximation based on Eq. ( ^d|) rather than on Eqs. ( ITbl ) and 
(ITTI) results from the fact, that Eqs. (P) and ( [lOl) treat the glass structure as it comes out 



from Eq. (H), thus avoiding the additional small-g approximation Cq ~ Cq=o- 

C. The stiff glass approximation 

In this section Eqs. (|l]) and (0) shall be considered for large coupling coefficients Vq^kp- Let 
us write symbolically Vq^kp = 0{l/ri), so that various quantities can be classified according 
to their power of the small parameter rj. For the function Cq, which enters Eq. (|^) for fq, one 
gets large values: Cq = Oil/rf). Thus 1 — fq = 0{ri). Stiff glass states are characterized by 
Debye- Waller factors close to unity, as is demonstrated by curve (a) in Fig. |1|. Therefore one 
concludes from Eqs. (1^), (|5c|) , ( |10b|) and (|10c|) that the renormalized coupling coefficients 
decrease towards zero in the limit rj ^ 0: 

Vqk, Wi'^ = 0{ii) , Vq,up, lyf = 0{ri') . (12) 

After eliminating the arrested glass structure as described in Sec.IIB, the remaining MCT 
equations (|4d| ) and (|) deal with a weak-coupling situation if the system is in a deep glass 
state. Figure |Tl| demonstrates how the coefficients Vr^^'^-* decrease with increasing packing 
fraction. Let us introduce for later use also the integrated coupling coefficients 

w^2= r w^'^^^dk . (13) 

Jo 

For the three packing fractions (p = 0.600 (0.567, 0.540) dealt with in Fig. ^ one finds wi 
= 0.39 (0.53, 0.64) and W2 = 0.14 (0.34, 0.78) respectively. Let us recall from Sec.IIB that 
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Vq^k and the derived quantities WJ: , wi specify the interactions of density fluctuations with 
the arrested amorphous structure, while Vq^kp, W^. ,W2 quantify two-mode-decay processes. 
One concludes furthermore that in the limit rj ^ the two-mode contributions to the 
kernel mq(u) get suppressed relative to the one-mode contributions, in particular: K^^\t) = 
0{r]), K^'^\t) = 0{rf). This explains the result shown in Fig. 12: for the packing fraction 



(/9 = 0.60 the one-mode contribution K'^'^'>" [uj] provides the dominant part of the total 
spectrum K"{uj). The role played by the two contributions K^^^ and K^'^^ is utterly different. 
The former yields the AOP and the latter provides a background spectrum. 

To explain the AOP a stiff-glass approximation shall be analysed which is indicated by 
superscripts (1). It is defined by dropping the two-mode contribution to rhq^uj) in Eq. (0). 
Substituting the result into the formula (p^^^cu) = [xi^H^) / Xq ~ 1]/^ one gets 

(l)^^\uj) = -l/[uj-nl/[uj + tiyq + nlm'^^\uj)]] , (14a) 

m(i)(cu) = ^«[0(i)(cu)] . (14b) 

In order to ease the discussion at the end of this section, a friction term proportional to 
Uq > has been included in the formulas. It is equivalent to complementing m^}^^"(uj) by 
a white noise spectrum. Unless emphasized otherwise, one may read the formulas with 
Uq = 0. Equations (0) define a special model for the MCT and all general theorems quoted 
in Sec.IIA apply. If one would treat the mode-coupling coefficients Vqk > in the Eq. (|5ED 
for the functional JF^^) as free parameters, Eqs. ( ^^ could exhibit glass transitions. For the 
discussion of such bifurcations the M by M matrix Cqk = {dJ-'q[f]/dfk} (1 — /fc)^, q,k = 
1,2,... M, which is called the stability matrix, plays an essential role. Glass transitions are 
characterized by the spectral radius E of matrix C to be unity, and for all other states one 
gets E < 1 p2|,^. One checks, that the stability matrix of the complete theory as defined 
by Eqs. (^) and (^) is the same as the one for the stiff-glass approximation, defined by 
Eqs. ([T^). Therefore Eqs. (|1^ with the coefficients Vqk defined in Eq. (|5BD do not exhibit 



glass-transition points anymore; in particular (j)^q\t -^ oo) = 0. Let us note also that matrix 
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C is equivalent to matrix V : Cgk = (1 — /g)V^fc/(l — fk)- Therefore also the spectral radius 
of V is E. Hence the resolvent S = {1 — V)^^ exists. 

The u = limit of Eqs. (|l^) yields the linear equation Y^ki^qk — ^fe)"^A,. ("^ = 0) = 
i J2k Vqk^k/^l- For the relevant case u^ = one concludes that the zero-frequency spectrum 
of the kernel vanishes. Therefore 

mW"(cu ^ 0) = Rqto^ , (15a) 

and this implies for the low frequency behavior of the correlator 

0«"(l^ ^ 0) = (7r/2) [6iuj - fif ) + 6iiu + fif )] + Rqcu^ . (15b) 

Hence, one finds Rayleigh's law: the scattering probability m^^"{Lu) of low-frequency 
phonons from static density fluctuations varies proportional to u"^. 

Let us simplify the stiff-glass approximation by embedding it into the generalized- 
hydrodynamics description. According to Sec.IHB this amounts to approximating m^^\uj) 
by its g = value, denoted K^^\uj). Similarly one should write for the so far not specified 
friction coefficient Ug = tCi'^. Let us denote the stiff-glass susceptibility resulting from Eq. 
( |14a|) within the generalized-hydrodynamics approximation by x^'^^(a;) = Xq{(j)''g\uj)-uj + l), 



so that 

Xf ^(^) = Xq^V [^' ~^l + ^^l [ir + K'^^u;)]] . (16) 

The problem is reduced to evaluating the kernel K^^^{uj) from the equation 

roc ^ ^ 

K^GA^^^ ^ _ / dkWt'^{l/ L - nil \u + nl{iT + K^'^^u)]]} . (17) 

JO L L J J 

Again, the general theorems of MCT apply. However, since the mode coupling coefficients 
have been altered by substituting Vq^ for Vqk, the stabihty matrix has been changed. There- 
fore it cannot be excluded that the spectral radius E reaches unity for the simplified theory 
. This would be an artifact of the generalized-hydrodynamics approximation and Eqs. (|TB[) 
and (|T7|) must not be used in such case. 
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The generalized-hydrodynamics approximation simphfies the Eq. ( |1 5a|) by the replace- 
ment of Rq by its g = limit to be denoted by Rq : K^'^^"{ijj ^ 0) = Rquj"^. The coefficient 
Rq is obtained by substituting Eq. ( |15b| ) into Eq. ( |10b|) . Using Eq. ( pE| ) one gets 



R, = p{S,lv^Yclh{l - hYl A-K^l - w,) 



(18) 



The Rayleigh spectrum is included in Fig. |1^. A huge magnification was necessary to make 
this contribution visible on the scale of K"{uj). 

One infers from Fig. ^ that for q > 2 the variations of the renormalized frequency Qq 
with changes of the wave vector q is suppressed relative to the ones of Qq. One can get an 
estimation of the k > qo contribution to the integral (0) by replacing Cl^ by some averaged 
value Cl. Figure suggests Q = 120 for ip = 0.60. Introducing u = rCl"^ and indicating the 
results of the specified estimation by a tilde, Eq. ([T6|) is equivalent to 



^q>qiiuj) ^ 4)iiu) = -1/ iU-ny iU + tl^ + ^^ K (uj) 



(19a) 



According to Fig. [TT] the k < qu contributions to the integral in Eq. (|T^ can be neglected, 
and thus this equation for the kernel simplifies to 



K 



SGA 



(cj) ^ K{uj) = Wi(j){uj) 



(19b) 



The solution of Eqs. ([T9|) for K and (j) reads 



kiu) 



uj,uj_ — ziio) + \ ziujV — Lot ■ \ ziujV — UJ 



2 _,.,2 



\2 _,.,2 



/(2c.) 



(20a) 



z{ijjY = ^{^ + iy)/^^ , 



u;± = 1 ± ^Jwl . 



(20b) 



The coupling constant Wi takes over the role of the spectral radius of the stability matrix. In 
order avoid artifacts of the various approximations leading to the results (|20|), the formulas 



can be applied only for wi < 1. Figure |T3| demonstrates that Eqs. ( pOf ) with z/ = describe 
the AOP of K"{uj) for the HSS with (f = 0.60 reasonably well. One could improve the 
description by trying better choices for Q, but this would not lead to any new insight. 
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Figure |TT] shows that also for the evaluation of the integral in Eq. ( |10c|) the contributions 
for k < q£, can be ignored. A leading estimation for the two-mode kernel can thus be 
obtained as K'^'^\t) = ti?20(i)^- Using Eq. (|19bD one gets therefore 

i^(2)(^) = {1/tt){w2/wI) I k{u - uj')k'\J)dJ . (21) 



The zero-frequency limit leads to a trivial integral with the result k'^'^'>"{ijj = 0) = r = 
{W2/ ^/wi)S/{?>tiVL). Going back to Eq. ( |10a| ) one concludes: the background term due to 
K^'^\t) could have been taken into account in its white noise approximation. This leads 
to the extension of the equation of motion by adding a friction term v = D? ■ t. Such an 
extension is obtained by including this term in the formulas, as was done already in Eq. 
( |14a|) and the following formulas. The dotted line in Fig. |13| demonstrates that thereby 



all qualitative features of K"{uj) are understood. A further improvement is obtained by 
dropping the white noise approximation for the correction term K^"^' . This is done by 
replacing cu + iu hj to + K^^^'^to) in Eq. (|20b|) . Figure |l^ demonstrates that thereby a more 



satisfactory treatment of K"{uj) is obtained. 

The periodic continued fraction for (piuj), which is defined by Eqs. (p^), can be related 
to a Hilbert-Stieltjes integral: j\ dx\J\ — x^ /{x — Q = 7r(— ^ + ^C, — IVC + !)• Thus one 
can express the normalized susceptibility xi'^) = [^0(^) + 1] i^ the form: 

X{uj) = / drp{i)x^{uj) , (22a) 

m = ^{u:l-m-ujl)/{27,w^) , (22b) 

Xt,{uj) = -f2V [tu^ - n"^^ + iuu] . (22c) 

For the stiff-glass states the AOP is obtained as a superposition of undamped-harmonic- 
oscillator spectra. The weight distribution p(^) for the oscillators with frequency y^^ 
extends from i7_ = uo^Vt to fi+ = uoj^VL. If the approximate description is extended so 
that two-mode interactions are incorporated as white-noise-background spectrum for the 
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fluctuating-force kernels, the results remain valid, but the oscillator dynamics has to include 
a Newtonian friction term, quantified by z/ > 0. A better description of the spectra for 
frequencies large compared to Q^ is obtained by acknowledging that the friction forces do 
not exhibit a white noise spectrum. This can be done by replacing iu in Eq. ( p2cD by the 
kernel K'^'^\uj) from Eq. (PTF). The glass instability for T ^ T^ or (^ ^ ^pc is connected with 



the approach of Wi to unity, i.e., with the threshold fi_ approaching zero. 

IV. SCHEMATIC MODEL DISCUSSIONS 

A. Models for anomalous oscillation peaks 

The simplest MCT models deal with a single correlator, say 0(t). The Eq. (^) remains 
valid with the subscript q dropped. Generalizing Eq. (^), the kernel m{t) is written as 
mode-coupling function, specified by a series with coefficients Vn > 0: 



m 



(t) = J' [0(t)] = y: vn<p{tr . (23) 



n=l 



The long-time limit / = 0(t -^ oo) obeys Eq. (|^): / = jF[/]/(l + jF[/]). For the glass states 
one can carry out the transformation discussed in Sec.IIB to get correlators (j){t) and kernels 
m(t) with vanishing long time limits. Equations (^ and (^) remain valid with the subscript 
q dropped. The new functional T is again a power series. For the renormalized Taylor 
coefficients one gets in analogy to Eqs. (BE) and (p^): {)„ = (! — /)""'"^[9"JF[/]/(9/'^]/n!, n = 



1,2,.... The limit of a stiff glass is obtained if at least one of the coefficients f„ becomes 
large. The number rj = 1 — f can be used as small parameter for the classification of terms. 

One finds {)„ = 0{ri"-). Let us change the notation to Wi^2 = "^^1,2 so that 

00 00 

w, = il- ff Y. ^^nP-' , W2 = il- ff E n{n - l)vr,r-' ■ (24) 

n=l n=2 

The stiff-glass approximation can be defined as in Sec.IIIC by approximating the kernel by 
the leading term m(t) = wi(j){t). Let us change the notation to = and Cl = Q, so that 

0(t) = / + (!- mt) , n' = fiV(i - /) • (25) 
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With K{t) = rh{t) one arrives at the pair of equations (pj|). 

Let us consider as an example the model specified by the functional JF[/] = vif + V2p. 
It was introduced as the simplest one which can reproduce all possible anomalous exponents 
of the general MCT p2| , [28[] . Liquid-glass transitions occur on a line in the Vi — V2 parameter 
plane. The line, where the long-time limit / jumps from zero to /'^ = 1 — A > 0, is a 
piece of a parabola with the representation vl = (2A — 1)/A^, V2 = 1/A^, 0.5 < A < 1. In 
Ref. [^ diagrams analogues to Figs. ^ and |^ can be found, which exhibit the evolution 
of glassy dynamics and the AOP upon shifting (fi,f2) from the weak-coupling regime to 
the strong-coupling one. As path in the parameter plane a straight line was chosen: f 1 2 = 
■^12(1 + f)) -^ = 0.7, e = ±1/4", n = 0, 1, . . . The full lines in Fig. |l^ reproduce three 
glass-state results. They are similar to what is demonstrated for the HSS in the lower panel 
of Fig. ^ The reason was explained in Sec.IIIC: these HSS spectra are described by the 
stiff-glass approximation and this approach yields the same formulas (0) as derived above 
for the schematic model. The dashed lines in Fig. |1^ represent the leading order description 
by Eqs. ( [T9| ) with u = and Wi,Q determined by Eqs. ( |2^ and (|25|). The dotted lines 



incorporate v = ri'^K^'^^"{uj = 0) evaluated from Eq. pT|). The results for the n = 4 spectra 
show that the extended stiff-glass approximation can describe the dynamics reasonably well 
for states which differ from the instability point by only 6.25%. It was pointed out |^ 



that the schematic-model solutions in Fig. |T4| fit qualitatively to the boson-peak scenario 
as reported for neutron-scattering data. It was shown ^^ that the AOP of the specified 
model can be used to fit the boson-peak spectrum of orthoterphenyl as it was measured 



by Raman scattering at T = 245i^ = T^ — 45K for frequencies between 100 GHz and 700 
GHz. 

Usually, a one-component schematic model is too restrictive to deal quantitatively with 
experimental data. But one can construct more elaborate schematic models with the inten- 
tion to mimic more features of the MCT. The perspective of such approach shall be indicated 
by results for the model used recently for the interpretation of scattering spectra of glassy 
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liquids [p9| , pl| -|3^ . The model extends the one of the preceding paragraph by introducing a 
second correlator, to be denoted by <f)''^{t). The equation of motion has the standard form of 
Eq. ( p!a]) with Qg, (j)g, and rriq replaced by Q^, (p^, and m^ respectively. The mode-coupling 
functional for m^ is characterized by a single coupling constant Vs > 0: 



m^(t) = t;,0(t)0^(t) . 



(26) 



The model was introduced for the description of tagged-particle motion in liquids 

The long-time limit /* = (f)^{t -^ oo) is obtained from Eq. @ as /'^ = 1 — l/{fvs)- The 
transformation to a new correlator and a new kernel with vanishing long time limits can be 
done as explained in Sec.IIB: Vt'"^ = fi^V(l - /"); 4>%t) = [4>'{t) - /"]/(! - f). From Eq. 



(|cD one gets 



m^it) = Ki'\t) + Ki'\t) 



(27a) 



Ki'\t)=u^it)+ws^'{t) 



Kf\t)=Vs4>{tW{t) 



(27b) 



u = {l-f)r/f 



Ws 



;i-r) 



;i-r)(i-/)// 



(27c) 



The stiff-glass approximation for the two component model requires (1 — /) and (1 — /*) 
to be small so that K^^\t) can be neglected compared to K^^^t). Let us denote the results 
by a tilde. The equations of motion, which specialize Eqs. (p!4D, read: 



(iu) = -1/ Lj- fi"7 Lu + ifiuj) + niK'iuj) 



(28a) 



K'iu) = w,(h'(u;) 



(p{uj) = fi M0(a;) . 



(28b) 



This result is equivalent to Eqs. (|T^) except that the friction constant iv is generalized to 
a friction function <f{uj). The solution can therefore be written in form of Eq. ( ^Oa] ) with 
z^uY and uj± replaced, respectively, by 



z'{uY = uj{uj + '^{uj))/n 



s2 



oo: 



1 ± ^A(^ 



(29) 
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In particular Eqs. p^) hold with the appropriate change in the notation. The susceptibility 
is a superposition with weight p'*(0 of the functions 

Xlioo) = -Q'y [uj^ - n'^^ + uj^iuj)] . (30) 

This is a harmonic-oscillator susceptibility where the interaction with the background is 
included via a friction function (p{uj), dealing with the AOP of the surrounding. 



Figure 15 shows solutions for three characteristic choices of Vf^. The results refer to the 



stiff glass-state discussed in Fig. for the label n = 0. The coupling f * = 10// was chosen 
so that 1 — f^ = 0.10. The stiff-glass approximation results are shown as dashed lines. 
For Q^ = O.lSfi the resonance of the second correlator is located so far below the AOP of 
the first correlator, that there appears an AOP of (^'^"{uj) quite similar as discussed for the 
one-component model. In this case (p{uj) in Eq. ( |30D only produces a renormalization of 
the frequencies ^^y/^ and a strongly suppressed background. For Q'^ = Q the resonances 
of Xf are shifted upward because of level repulsion and the hybridization yields a broad 
background extending from the threshold of the AOP of the first correlator to the AOP 
position of the second one. For i7* = 0.6i7 the spectrum (p'^'^uj) exhibits an AOP whose 
low frequency threshold and maximum position are close to the ones for the AOP of (f)"{uj). 
The hybridization causes a suppression of the high frequency spectrum. Therefore the AOP 
in (f)''^"{uj) is more asymmetrical than the peak in (p'^uj). The hybridization results have 
similarities to the ones discussed in connection with Fig. ^. Indeed, it was explained in 
Sec.IIIC that the phonon modes are influenced by the large- wave-vector modes which built 
the AOP, but that there is no feedback of the phonons for wave vector q < qo on the AOP. 
This is the same situation as treated by the specified schematic two-component model. 

B. Random oscillator models 



To get more insight into Eqs. (p2D, the following problem shall be considered: evaluate 
the averaged dynamical susceptibility x{^) of ^in ensemble of independent harmonic oscilla- 
tors. The oscillators are specified by their mass /i and by their frequencies fi(0 > 0. The 
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latter depend on a random variable $,. Its distribution shall be denoted by p{C); p(0 ^ 
0, /p(0'^^ = 1- It is no restriction of generality to assume ^^(0^ = a + P^, /^ > 0; but C, has 
to be restricted from below to insure stability. Let us choose as minimum for C, the value 
-1. Thus one can write fi(0^ = ^o[l + ^ + 2v^^]) so that ^o > defines the frequency 
scale and w,0 < w < 1, characterizes the minimum frequency Q- = r2o(l — V^)- Denoting 



averages by bars, A{^) = J A{^)p{^)d^, the quantity of interest is 



Xiuj) = -pI [^2 _ i](^)2] . (31) 



Let us define a characteristic frequency fi > by fi^^ = Vt{^)^'^. It specifies the static 
susceptibility x = p/Q"^. If brackets denote canonical averaging, defined with respect to 
the oscillator Hamiltonian H^ = [{P^/fi) + fiQ{^)'^Q'^]/2, one gets for the fluctuations of 
the momentum (P^) = /i^f ^ and of the displacement {Q"^) = v'^/Q{^)^, where v denotes the 
thermal velocity. The time evolution of some variable A = A{Q, P) can be written as usual in 
terms of a Liouvillian C : A{t) = exp{ia)A, where iCA = {dA/dQ)P/fi-{dA/dP)fin'^{OQ- 
To embed the problem into the standard framework of correlation-function theory a 



scalar product shall be introduced in the space of variables A, B . . . by {A \ B) = {A*B). 
The vectors | Q) and | P) are orthogonal and the normalization reads {Q \ Q) = v^/Q'^, {P \ 
P) = n'^v'^. The Liouvillian is hermitian. The displacement correlator shall be defined by 

0(t) = {Q{t) I Q)/iQ I Q) . (32) 

Its Fourier-Laplace transform can be written as Liouvillian-resolvent matrix element: 
(f){uj) = {Q \ [C — uj]^^ I Q)/{Q I Q). This quantity can now be represented within the 
Zwanzig-Mori formalism as a double fraction [35|: (f){uj) = — 1/[^ — ^"^/[uj + Q'^m{uj)]]. The 



memory kernel m{uj) is the Fourier-Laplace transform of the fiuctuating-force correlator 

m{t) = (F(t) I F)/{v^n^fi^) . (33a) 

Here F is the projection of the force dtP = —piQ'^{C,)Q perpendicular to | Q) and | P): 

F = fi[Q^-Q\0]Q . (33b) 
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The time evolution in Eq. ( P3a|) is generated by the reduced Liouvilhan £', F(t) = 
ex-p{iC't)F, where C = VCV and V denoting the projector perpendicular to | Q) and 
\P). The susceptibility is connected with the correlator as usual, x(^)/x = (^0('^) + 1); 

X{uj) = -At/ [uJ^ - fi^ + ujQ^m{uj)] . (33c) 

This exact representation of x{^) i^i terms of kernel m{uj) is the analogue of Eq. (0). 
The essential point in the MCT is the approximation of the kernel m{t) as mode- 



coupling functional. The procedure |^, consists of two steps. Firstly, one reduces 
F to the projection on the simplest modes contributing, and these are the pair modes 
I Q • ■ ^(J^) ^1 Q{'t)OiQ^ I QO~^iQ^ I ^)- Secondly, one factorizes averages of products 
into products of averages: {Q{t)^ \ Q^) — > {Q{t) \ Q)C,'^. As a result one finds 
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m{t) = wi<P{t) , w, = i m^in\0) - 1] le ■ (34) 

These formulas are the analogue of Eqs. (^) and (||). They allow the approximate evaluation 
of x{'^) from the given input information Vl? and Wi. 

The Eqs. ( |33c| ) and ( |5^ are equivalent to Eqs. ( [TP| ) with Vt = Vt = Vt^ and therefore the 



approximation for x{'^) can be written as noted in Eqs. (^). Thus, the presented MCT 
delivers an approximation for the oscillator susceptitiblity (J3l|) in the sense that the general 



distribution p{S,) is approximated by p(^) = 2-^/1 — E,"^ /n. If p{i) = 2-^/1 — .^^/tt is chosen, 
one can check that fi = i7o and w = wi. In this case MCT reproduces the exact result. Let 



us add, that the naive factorization for the kernel, (-F(t) | F) ~ iJ?{Q{t) \ Q)[Vt'^ — Vt{C)'^Y^ 
would not reproduce the exact result for the specified example; rather one would obtain Eq. 
(P^) with an overestimated wi = w + w"^. 

V. CONCLUSIONS 

Within mode-coupling theory (MCT) a critical temperature Tc and a corresponding 
critical packing fraction ip^. were introduced as the equilibrium-thermodynamics parameters 
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characterizing the evolution of glassy dynamics. For silica, for example, Tc is near 3300K 
P7| , and therefore all experiments quoted for this system in Sec. I deal with T <^Tc states. 
In this paper primarily states are studied where T is so far below Tc and (p so far above (pc 
that structural-relaxation phenomena do not dominate the dynamics within the window of 
interest. These states are referred to as stiff-glass states. The dynamical window considered 
is the one of normal-condensed-matter physics, i.e. spectra are discussed within a two- 
decade regime for the frequency uo around and below the Debye frequency uj£,. For the 
stiff-glass states the spectra of the a-relaxation process are located at frequencies smaller 
than ujd/^^^ and therefore it does not matter whether or not the quasielastic a-peaks 
of the spectra are treated as elastic ones. Hence it is legitimate to use the basic version 
of the MCT which treats the crossover near T^ as a sharp transition to an ideal glass at 



Tc. The derivation of the MCT formulas, in particular the one of Eq. pB]) for the mode- 
coupling coefficients, is based on canonical-averaging properties. For temperatures below the 
calorimetric glass-transition temperature T^, the system is in a quenched non-equilibrium 
state. From a rigorous point of view the application of MCT is therefore restricted to the 
regime T > Tg. However, experiments on high-frequency sound and on the boson peak do 
not indicate anomalies for T near Tg. Thus it seems plausible that the results of the present 
paper can be used also for an interpretation oi T < Tg data. 

A major finding of this paper is that there are "peaks" of the density-fiuctuation spectra 
for wave vectors q exceeding about half of the Debye-vector qd, which are quite different 
from what one would expect for phonon resonances in liquids or crystals. These peaks, 
which we refer to as the anomalous-oscillation peaks (AOP), show the properties of the 
so-called boson-peak listed in Sec. I. First, as shown in Fig. 7, the position up of the peak 
maximum is several times smaller than un. In this sense the AOP is due to soft modes. 
Second, according to Eqs. (p2D, the AOP is a superposition of harmonic-oscillator spectra, 
where different oscillators are specified by different frequencies. This can be shown since 
there is a well-defined strong-coupling limit of the MCT equations, referred to as stiff-glass 
limit, where the equations of motion simplify so much that all features of the AOP can be 
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worked out by analytical calculations (Sec.IIIC). In this limit the continuous spectra are 
purely inhomogeneous ones. Third, there is a lower cutoff Q^ of the frequency distribution. 
This causes the low frequency wing of the AOP to decrease more steeply with decreasing 
uj than expected for a Lorentzian. The high frequency wing extends further out than the 
low frequency one so that the AOP is skewed. Fourth, as the packing fraction decreases 
also the frequencies up and fl- decrease, and simultaneously the intensity of the spectrum 
increases. This is shown in the lower panel of Fig. ^ and in Fig. |10[ The critical point is 



characterized by the threshold Q^ approaching zero. In this sense one concludes, that the 
evolution of the AOP is related to the dynamical instability predicted by MCT for ip = ipc 
or T = Tc. In order to understand a further property of the AOP one has to acknowledge 
that the leading correction to the stiff-glass results introduces a damping for the oscillators. 
It is due to the decay of an oscillator mode into two modes caused by anharmonicities. A 
simplified treatment of this phenomenon only leads to a modification of the formulas by 
the introduction of a friction constant u in Eq. ( |22c| ). Thus the peak is superimposed on a 
flat background and the sharp thresholds are changed to some smooth but rapid crossover. 
Not much is modified in the center of the AOP spectrum provided u is not too big. But 
with decreasing ip or increasing T, the ratio z//fi_ becomes much larger than unity so that 
the oscillators of low frequency get overdamped. As a result one obtains the explanation 
of the fifth property, namely a central relaxation peak is formed if {ip — ipc)/Vc is about 
10%. In this case the AOP merely appears as a shoulder of the quasielastic spectrum as 
shown by the n = 4 curve in the lower panel of Fig. |^. Shifting the parameters even closer 
to the instability point, the AOP gets buried under the wing of the quasielastic peak. The 
elementary formulas (p!9D and (^) describe this feature of the MCT solutions reasonably well 



for Lp approaching ipc up to about 5%. The cited results are quite general and are obtained 
even for the simplest schematic MCT models, as is demonstrated in Fig. |1^. Sixth, in a 
further refinement of the description one acknowledges that the glass compressibility has a 
wave-vector dependence. This enters in the form of a characteristic frequency Vtg, which can 
be considered as a bare phonon dispersion of the amourphous solid. It exhibits a maximum 
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for q near qa and a minimum near the structure-factor-peak position q^iax which in turn is 
near 2g/). This lig-versus-g curve is similar to the one for the characteristic frequency Vtq 
of the hquid which plays an essential role in the MCT equations (p]a| ) and (jT^). However, 
the oscillations of the Debye- Waller factor, shown in Fig. [^, imply via the renormalization 
formula (HH), that the ratio of the maximum to the minimum frequency is much smaller 



for Vtq than for Vtq. This can be inferred in detail by comparing the full line in Fig. ^ with 
the dashed one. Therefore, the maximum position of the AOP is only weakly g-dependent. 
Summarizing, we suggest that the MCT of the AOP provides the basis of a first-principle 
explanation of the so-called boson peak. 

The formulated theory for the AOP has a transparent interpretation. MCT explains in 
the first place the formation of an effectively arrested density distribution. This amorphous 
structure is characterized by the same quantity used to characterize crystalline structures, 
viz. by the Debye- Waller factor fq. Within the same formalism which leads to /g, the 
equations of motion for the density fluctuations of this structure are obtained (Sec.IIB). 
Such unified treatment of the glass structure and its dynamics should be a feature of every 
microscopic theory since it is the same array of particles which forms the frozen structure 
and which carries the fluctuations. This unified treatment is especially important if one 
intends to study the dynamics near the instability limit of the structure. The equations of 
motion describe the decay of fluctuations into pairs and also the scattering of fluctuations 
from the arrested structure. For the stiff-glass states one finds the latter processes to 
overwhelm the former, as discussed in connection with Eq. (p!2D. In the stiff-glass limit one 
finds that the particles are localized in their cages and harmonic oscillations of the particles 
with their cages are a good description of the relevant modes. In this extreme limit the 
total susceptibility is the one of a distribution of independent-oscillator responses, where 
the distribution of oscillator frequencies is caused by the distribution of sizes and shapes of 
the cages. MCT provides an approximation of the distribution of the frequency squares, Eq. 
and characterizes the drift of the distribution with changes of control parameters. It 



provides also results for the corrections to this limiting result, namely the appearance of 
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homogeneous line broadening due to mode decay and coupling of the oscillations leading to 
weak wave-vector dependencies of the AOP parameters. 

One achievement of MCT is the possibility to explain homogeneous line-broadening 
effects via a golden-rule mechanism as is suggested by Eqs. (Q) and (|5^). This aspect was 
used above in connection with the evaluation of the oscillator damping v due to two-mode 
decay and also in connection with the derivation of Rayleigh's law, Eq. ( [15a| ). However, 
the interpretation of Eqs. (|5d|) and (^ in the spirit of a golden rule is quite misleading, if 
the transition probabilities V are so large that the self-consistent solutions are qualitatively 
different from the ones obtained by a lowest-order approximation. In such cases the formulas 
can lead to an approximation theory for an inhomogeneous "line width" phenomenon. It 
was shown explicitly in Sec.IVB that MCT provides an approximation approach towards 
this phenomenon, and there is an example for which MCT reproduces the exact result for 
the inhomogeneous spectrum. The MCT for the AOP is based on the fact that this theory 
can handle homogeneous and inhomogeneous spectra within the same framework. 

A side remark might be helpful. According to MCT all structural relaxation features 
are independent of the details of the microscopic equations of motion p5| , ^ . Therefore 
the existence of an AOP is of no relevance for understanding glassy dynamics or the glass 
transition. But the AOP provides an interesting piece of information on the arrested glass 
structure. The AOP is the result of a mapping of the cage distribution on the frequency axis. 
The mapping is done from the configuration space on the time axis via Newton's equations 
of motion and canonical averaging followed by a Fourier-cosine transform to get a spectrum. 

The arrest of density fluctuations at the ideal liquid-glass transition is driven by the 
ones with a wave number q near the structure-factor-peak position gmax ~ 2g£), since for 
these wave numbers the liquid compressibility xJ oc Sq is large. The compressibility of 
simple dense liquids for q < q^ is very small and therefore excitations with wave vectors 
from this domain, which corresponds to the first Brillouin zone of the crystalline phase, are 
not important for the evolution of structural relaxation and the glass transition. For the 
same reason one concludes that scattering processes of density fluctuations with q < qo are 
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irrelevant for the formation of the AOP, as is demonstrated in Fig. 0. The soft complexes 
which cause the AOP are constructed from fluctuations with wave numbers near and above 
?max- These conclusions are based on the MCT results for the Debye- Waller factors fg. In 
order to produce the spectral peak in Fig. |T2], fg has to be that large as shown by the 
uppermost curve in Fig. ^ This curve corresponds to e = {ip — ipc)/fc = 0.16. In Fig. ^ of 
Ref. I^Sf a measurement of fg for a hard-sphere colloid is documented for 3 < qd < 13 and 



e = 0.11. Since these experimental findings are close to the curve (a) in Fig. |l]we argue that 
the MCT results on the glass structure are in reasonable accord with the experimental facts. 
Let us emphasize, that the above reasoning refers to densely-packed systems of spherical 
particles. Obviously, in more complicated systems, like silica, the cages are not so tight 
like in a HSS. Therefore one can expect the soft configurations to be more subtle than 
discussed in this paper. The results in Sec.IVA for the two-component schematic model 
indicate, that the AOP can be more structured than obtained for the HSS. For systems with 
a low coordination number one can also expect intermediate range-order effects to play an 
important role. They enter, e.g., the coupling vertices in Eq. (PE) via the prepeaks of Sq 



\ . Whether MCT can handle the microscopic features leading to the AOP in complicated 
systems such as ZnCl2 is unclear at present. In particular, it is unclear whether MCT can 
contribute to understanding why the boson peak is more pronounced in strong glass formers 
like silica than in fragile glass formers like orthoterphenyl. 

The preceding interpretation of the AOP suggests that these peaks appear in the spectra 
of all probing variables coupling to density fluctuations of short wave length. But different 
probing variables will weight the oscillating complexes differently and therefore the shape 
of the AOP and the position up of the peak maximum will depend somewhat on the probe. 
Let us consider two examples. The first one deals with the tagged particle correlator ^^(t) = 
(p|(t)p|-*). Here p|(t) = exp(zgF(t)) denotes the density fluctuation of a marked particle with 
position vector r(t). The spectrum (pg'iuj) determines the incoherent neutron-scattering 
cross section. An exact equation for this quantity has the same form as Eq. (|T^ with (pg, mq 
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and fig replaced by 0^, m* and ^2^^ = t>^g^, respectively. The essential MCT equation is again 
the representation of the kernel as mode-coupling functional m^(t) = Y,kp^q,kp4'k{t)4'p{i) ■ 
The coupling coefficients V^^^^ are determined by the structure factor Sq. The mentioned 
equations have been derived and solved for the Lamb-Mossbauer factors /* = 0g(t ^ oo) 
of the HSS in Ref. |^. Details of the discretization can be found in Ref. |^T|. We have 



solved the cited equations for 0q(t), and Fig. |l^ exhibits fluctuation spectra for the packing 
fraction ip = 0.60 for three wave vectors around the structure-factor-peak position. The 
shape of the peaks is only weakly g-dependent and the intensity varies nearly proportional 
to g^. This finding is in agreement with neutron-scattering results EOl and with molecular- 



dynamics-simulation results for ZnCl2 ||T6[ and silica |50l . 

The second example is related to the velocity correlator "^(t) = (r(t)f)/(3f^). For a 
harmonic system its spectrum determines the density of states g{uj) = 2'$"{uj)/tt, normal- 
ized by J^ g{uj)duj = 1. The velocity correlator can be extracted from 0^ |^. If one 
introduces the kernel m^'^\t) = liiag^Q q'^m^ (t) one gets '^{to) = —l/[uj + v'^m^^\uj)]. The 
localization length r^, defined via the long time limit of the mean squared displacement 
{6r'^{t -^ oo)) = 6r^, determines the small-g asymptote of /^ and the pole of the kernel 
in the glass: /^ = 1 — (gr^)^ + 0{q'^), m^^^uj) = — l/(u;r^) + 0(u;°). One finds that the 
density of states vanishes proportional to u"^ for small frequencies as expected for an elastic 
continuum g{uj) = go ■ cu"^ + 0{uj^), g^ = m^'^^"{uj = 0)(2rj/7iv'^) [^]. Figure exhibits a 
result for the HSS. The tu^-law is obtained only for frequencies below the threshold Q- of 
the AOP. For larger frequencies the density of states is enhanced relative to the asymptotic 
law. The enhancement reaches a maximum of a factor of about 5 near the position of the 
AOP maximum, and it is about a factor 3 at the maximum of g{uj). For larger frequencies 
g{uj) gets suppressed as required by the normalization condition. The found enhancement 
phenomenon is in qualitative agreement with the experimental results reported for silica in 
Ref. |]T5[ and with the simulation data in Ref. [|T^. There is a consistency problem for the 



MCT. The prefactor in the tu^-law should be g^ = uj]J^[l + 2{ujd/^'£)Y], where u;^ denotes 
the Debye-frequency for transversal sound. Since go was calculated without any explicit 
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reference to transversal sound, the approximations underlying MCT will lead to go 7^ g'o- 
We did not study the solutions of the MCT equations for transversal excitations |2yj2^ in 
order to calculate uj'^. Therefore we do not know the size of the error go — g'o- But if one 
estimates ujd = ^'d 0116 gets g'^ = 10.2 10~®, which is close to the value go = 9.2 10"*^. This 
suggests that transversal excitations are taken into account to some extend implicitly in the 
formulas for m^Jt). 

It was already suggested earlier in connection with a discussion of soft-configuration 
models for glasses that the boson peak should be understood as result of quasi-harmonic 
oscillations of the system characterized by some distribution of oscillator potentials |5l|-|54 



Obviously, the present theory is consistent with these phenomenological approaches. In Ref. 
Q Raman spectra of glassy systems have been interpreted as a superposition of oscillator 



susceptibilities analogous to what is formulated in Eqs. (]22|). But there are two qualitative 



differences between this fit procedure and the present theory. In Ref. |^^ the distribution 
p{C,) is taken as temperature independent, while Eq. ( p2b| ) for p(^) describes the softening 
of the glass structure upon heating and, in particular, its instability for T reaching T^. In 



Ref. [^ the T dependence of the spectra is introduced by replacing the damping constant 
z/ by a Debye-function quantified by a temperature-dependent relaxation time. This visco- 
elastic theory leads to a Debye peak as quasi-elastic spectrum. Equations ( p2|) do not lead 
to a quasi-elastic Debye spectrum as was explained in connection with Figs. |^ and |^. In 
Refs. p5| , p6| the effective-medium theory for percolation problems is modified to a theory 
for the displacement susceptibilities of a disordered harmonic lattice. For the susceptibility 
an expression similar to Eq. ([T6|) is obtained where kernel K^'-^^^u) also describes the self- 
consistent treatment of phonon scattering by the disorder. Even though the equation for 
K^'-^^{ij) in Refs. |^5|,0 is quite different from Eq. (p!7D, the solution looks similar to the 



dashed line in Fig. 13. However, the Rayleigh contribution in Ref. [56| is about 10 times 



larger than the result based on Eq. (|T8|). It was criticized |^ that in Refs. |55| , p6| the 
boson peak is constructed from fluctuations with wave vectors q < Qd since thereby the role 
of long-wave-length fluctuations is overestimated. Hence the cited harmonic-lattice theory 
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does not appear to be compatible with the theory studied in this paper. 

Sound is obtained due to the interplay of inertia effects and stresses, which are built 
up due to compressions. The interplay is governed by the conservation laws for mass and 
momentum. The low-lying sound excitations interfere with other low-lying modes like 
structural relaxation and the oscillations building the AOP. Sound modes with wave vector 
q < qo and their interactions are not important for the explanations of structural relaxation 
pO| . In connection with Fig. |ll] it was shown, that they are neither relevant for the 



explanation of the AOP. Therefore sound can be discussed within the standard procedure of 
acoustics by introducing a modulus. This procedure comes out within MCT as was explained 
in connection with Eq. (|^), where [1 — ujK{ujy\ is proportional to the complex modulus. 
Consequently, all results on high-frequency sound, discussed in this paper, are implications 
of the preceding results for K{uj). In this sense one concludes that high-frequency sound is a 
manifestation of the AOP. Let us contemplate the scales for frequency uj and wave number q 
in order to be able to correlate the MCT results with some data. The boson-peak maximum 
observed for silica at 1 THz [Q shall be compared with the maximum position uop = 75 for 
the (fi = 0.60 results for the AOP. The best resolution width Fexp achieved by the recent 
X-ray-scattering experiments is 1.5 meV, i.e. in the units of this paper Fexp ~ 30. The 
structure factor peak position gmax = 1-5 A~^ for silica is to be compared with the value near 
7 for the HSS. Hence the wave vector unit used is about 2nin~'^. Scattering experiments 
with the resolution Fexp have been done for silica for q between lnm~^ and 4nm~^ as can 
be inferred from Ref. |^ and the papers quoted there. X-ray-scattering experiments with 
larger q are done with a resolution considerably worse than the cited Fgxp- Thus the following 
discussion shall be restricted to wave vectors between q = 0.5 and q = 2 ^ Q'd/2. This wave- 
vector interval corresponds to the interval for the sound frequency Clg = v^oq between about 
40 and 150 as is shown in Fig. |^. 



The first property of high-frequency sound follows from the upper panel of Fig. |T0. 
On the resolution scale Fexp the sound-dispersion law is Q^^^ = Qg = v^oq- Here the 
sound speed foo is the one determined by the g — > limit of the glass compressibility Xq. 
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The frequency-dependence of the reactive part of the modulus l — uK'^u) imphes deviations 
from the strict hnear law for Q^^^. According to Fig. |^ the deviations are predicted to occur 
on a 10% level. Thus they should be measurable if the resolution Fexp could be reduced by, 
say, a factor 5. Here some reservation has to be made. MCT does not contribute to the 
understanding of the structure factor Sg, rather it takes this quantity from other theories. 
Errors in Sq will cause errors in the MCT results. It is notoriously difficult to calculate the 
small-g behavior of Sq. The Verlet- Weiss theory yields different results for Sq than the used 
Percus-Yevick theory; and this causes also a small-g behavior of fq which differs from the 
one exhibited in Fig. [l| |^. Therefore the Verlet-Weiss theory will also lead to different Qq 
which might change the results for Q^'^^. It is not known how reliably the Percus-Yevick or 
the Verlet-Weiss theory describe the structure for packing fractions as large as 0.60. 

If one smears out the spectrum K"{uj) shown in Fig. 10 with a resolution curve of 
width Fexp, one gets a result which is nearly cu-independent within the dynamical window 
of interest: v'^K"{ijj) = 7. This explains the second property of high-frequency sound 
reported for the X-ray-scattering results of silica [^,0 and other systems Qj^j^l : the sound 
damping exhibits the hydrodynamic wave vector dependence Tq = ■yq'^. Consequently, one 
can describe the whole measured spectrum by the damped-oscillator formulas ( pa] ) and ( pb[ ) 
imj^J^, albeit up to some background. The latter appears as white if viewed with resolution 
Fexp- Naturally, it is difficult to separate this background from the one caused by other 
effects of the experimental setup. If one acknowledges the frequency dependence of K^u), 
exhibited in Fig. [10], one concludes, that the formula F^ oc g^ is oversimplified. Figure ^ 
demonstrates the prediction, that a reduction of Fexp by a factor of three should be sufficient 
to detect an increase of Tq above the F^ oc ^2™^^^ asymptote if q varies between 0.8 and 2. 

A crucial experimental finding is that the damping parameter Tq does not vary much with 
changes of temperatures. This shows [^,0 , that the sound damping mechanism cannot be due 
to anharmonicity-induced mode decay as known for phonons in crystals nor due to structural 
relaxation effects as studied in Brillouin-scattering spectroscopy of glassy liquids. The fast 
sound detected by neutron scattering [^ and molecular-dynamics simulation [^ in water. 
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for example, occurs in a dynamical window where water exhibits its a-relaxation process but 
no vibrations underlying a boson peak p0| , |61| . Therefore the fast-sound- damping in water 



depends appreciably on temperature and this dependence can be described reasonably within 
a visco-elastic model [^ . The insensitivity of the high-frequecy-sound damping on changes 



of control-parameters like temperature or density is indeed the third specification of the 
MCT results as shown in Fig. ^ In agreement with the assumptions of the phenomenological 



theories in Ref. [p2| , |54| MCT explains the damping to be due to absorption of the sound 
mode by the oscillations building the AOP. The AOP depends on control parameters, as 
explained in connection with Figs. and |l^. Therefore F^ is not strictly independent of 
control parameters. However, changes of density or temperature primarily redistribute the 
spectrum of K" and thus the spectrum in the centre of the AOP does not change much. 
But, if the resolution Fgxp could be reduced, a more subtle prediction could be tested. For 
j^max j^gg^]^ ]^4Q^ Yq is (/^-independent even if the system is driven as close to the critical 
point as shown by curve n = 4 in Fig. |T0|. For smaller f]™^^ the width F^ increases with 



decreasing ip, and this is due to the appearance of the quasi-elastic relaxation peak of K"{uj). 
But for larger f2™^^, the damping constant should decrease with decreasing yj; in this case 
the softening of the system reduces the density of states for high-frequency-sound-decay 
processes. Let us emphasize that all MCT results discussed in this paper are solely based 
on the wave-vector and control-parameter dependence on the structure factor 5*^. The 
explanation of the modulus and its drift with control-parameter changes is specified semi 
quantitatively by the three numbers only, which are specified in connection with Eqs. (^OD 
and (pTI). Therefore our results are predicted to be valid for all systems with a structure 
factor similar to that of our HSS model, in particular for Lennard- Jones systems or van der 
Waals liquids like, e.g., orthoterphenyl. 

The AOP oi K"{ijj) causes via Kramers-Kronig relations a frequency dependent reactive 
part K'{ijj). This implies that the sound resonance cannot exhaust the spectrum. Contrary 
to what holds for hydrodynamic sound, there must be a background spectrum. This is the 
fourth feature specified for high-frequency sound in Sec. I. Obviously, the detection of such 
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background is difficult in view of many other reasons producing backgrounds for the experi- 
mental scattering signals. However, as explained in connection with the dashed lines in Figs. 
|T3| and |14|, MCT implies a fifth property: for T <^ Tc or ip ^ ipc, there is an effective low- 
frequency threshold Q- for the background. Such threshold can be used to discriminate the 
background due to density-fluctuation dynamics from the one due to experimental artifacts. 
There is a mathematically equivalent manner to formulate the physics of the background, 
which is better adopted to the present problem [|l^,0. The dynamical structure factor 
S{q,u!) = Sq(f)g{uj), considered as function of q for fixed frequency u, represents the average 
of the square of the density-fiuct nation Fourier components which oscillate with frequency 
u. The coherent contribution to these fiuctuations leads to a peak at q^ = uj/voo- There 
is only a small contribution for q < q^j, since it is very difficult to excite long- wave-length 
fiuctuations in densely packed systems. But there is a structureless background for q > q^j 
extending to high values of q. It is caused by the large-g density fiuctuations produced by 
the distortions of the wave front due to the arrested amorphous glass structure. Figure p!8| 
exhibits a MCT result for the HSS, which is in qualitative agreement with the simulation 



results reported in Ref. [ITU 



All specified MCT results for high-frequency sound and the AOP can be described well by 
the combination of the elementary formulas (|^) and (|20| ) with K{uj) ~ K{uj). If regimes are 
considered, where relaxation can be ignored completely, one can use Eq. (pOb|) with u = 0. 
In this case, only the two parameters Cl and Wi need to be specified in order to quantify 
the result. Function K{uj) replaces the damping parameter 7 of the hydrodynamics-theory 
result, Eqs. (J^). Introducing the third parameter z/, the range of applicability of the results 
can be extended so, that structural-relaxation precursors are included. We suggest to use the 
cited formulas for an analysis of inelastic-X-ray-scattering data for high-frequency sound 
and of data for the evolution of the boson peak in glassy systems. 

The derivation of the MCT equations is definitive and leads to a well-defined model for 
a non-linear dynamics. The point of view adopted in this paper is the following: results 
for the model are derived to provide explanations of previously unexplained features of the 
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dynamics of liquids and glasses and to predict new results to be tested by future exper- 
iments. However, the "approximations" leading to the mode-coupling expression for the 
fluctuating-force kernel, Eqs. (|^), are uncontrolled and therefore the range of validity of 
MCT is not understood. Let us conclude this paper by pointing out four open questions 
concerning the foundation of MCT, which are of particular relevance for the study of vibra- 
tional excitations. The first problem concerns the absence of any influence of transversal 
excitations on the dynamical structure factor. The results (|15a|) and ( |T8|) for Rayleigh's law 



account for the scattering of a longitudinal phonon into some other longitudinal wave, while 
the expected contribution due to conversion into transversal sound waves is missing. Simi- 
larly, the spectrum of the AOP is due to longitudinal excitations only, while Horbach et al. 
p!8| report that transversal modes influence the boson-peak. The second problem concerns 



the mode-coupling approximation for the fiuctuating-force kernel mq{t) for short times. For 
a Lennard- Jones system Eqs. (pa]) and (pBD yield an overestimation of r7ig(t = 0) compared 
to the values known from Monte-Carlo results. A procedure for eliminating this problem 
was suggested in Ref . p2[ , but it is unclear whether it can be used for supercooled systems. 



One also expects that the mode-coupling kernel should be complemented by some regular 
term [^. In a simple treatment this would lead to an additional friction term i'g(f)g{t) in Eq. 



(p!a|). For a dilute system Uq could be calculated in the binary-collision approximation. The 
third unsolved problem is the evaluation of Ug for the dense systems under consideration. 
Such friction term would imply an additive correction to the friction constant u in Eqs. 
(pop and (0). If the resulting u would be too big, the AOP could disappear in favour of an 
quasi-elastic spectrum due to overdamped oscillations. Whether this happens might depend 
on the structure of the system. As fourth unsettled question a cut-off problem has to be 
mentioned. For the large packing fraction ip = 0.60, the results of our calculations would 
change somewhat if the used cutoff wave vector q*d = 40 was increased. This is due to the 
slow decrease towards zero of the direct correlation function Cq of the HSS for q tending to 
infinity. Since the introduction of a cutoff is equivalent to a softening of the hard-sphere 
potential it is conceivable that in the stiff-glass limit our model is rather a model for argon 
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than for hard spheres. It would be desirable to examine whether this cutoff dependence 
disappears if a conventional regular interaction potential is used. A serious bottleneck for 
such examination is the necessity to obtain reliable results for Cq for strongly supercooled 
liquids. 
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FIGURES 
FIG. 1. Debye- Waller factor fg (crosses) and one tenth of the Percus-Yevick static structure 

factor Sq (lines) of a hard-sphere system (HSS) at the packing fractions (p = 0.600 (a), 99 = 0.540 

(b), and the critical packing fraction ip = 0.516 (c). The inverse of the particle diameter, 1/d, is 

chosen as unit for the wave vector q. 

FIG. 2. Density correlation functions ^q{t) of a HSS as a function of time t for the wave numbers 
q = 3.4 and q = 7.0. The curves refer to the packing fractions tp = 0.6 and tp = (/^^(l it 10~"'^) 
with n given in the figure. Here ipc ~ 0.516 denotes the critical packing fraction. The curves with 
label c are the solutions at the critical point, which approach the long-time limits 734 = 0.36 and 
/7.0 ~ 0.85. The units of length and time have been chosen here and in all the following figures so 
that the hard-sphere diameter d = 1 and the thermal velocity v = 2.5. 

FIG. 3. Fluctuation spectra ^g(w) of the correlation functions shown in Fig. |2[ 

FIG. 4. Some of the correlators from Fig. |2| on a linear time axis. The full lines refer to 
glass states (n = 3 : ip = 0.567; n = 4 : 0.540) with the arrows indicating the long-time limits 
fq = liuit^Q^q{t). The dashed lines with label c exhibit the critical decay {(p = 0.516) and the 
lowest dashed curves refer to the liquid state for n = 4 ((/? = 0.549). 

FIG. 5. Fluctuation spectra ^q{uj) of the correlators shown in Fig. H. 

FIG. 6. Spectra ^^(f-i-') (solid lines) of a HSS at packing fraction ip = 0.60 for some wave 
numbers q. The dashed lines show the generalized-hydrodynamics approximation. For q = 1.0 
and q = 3.4 the dotted lines show the generalized-hydrodynamics approximation with fq replaced 
by fq=0- 
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FIG. 7. Frequency Qg = vq/y/^ (dashed line), renormalized frequency Qg = ^g/-\/l — fq (full 
line) and position of the global maximum 0,^^^ of the spectrum (diamonds) of a HSS at packing 
fraction ip = 0.60 as a function of the wave number q. For wave numbers at which two separate 
maxima of the spectrum can be identified the frequency position of the peak with the lower intensity 
is marked by the open circles. The vertical bars mark the frequency intervals where ^l'(w) exceeds 
half of the maximum intensity of the spectrum. The arrows point to the positions of the Debye 
wave number qD = (SGvry?)^''^ = 4.08 and the Debye frequency uid = VooQd = 309 corresponding to 
the high-frequency sound speed Voo = 75.8. For q = 3.4 (q = 7.0) one gets ^3.4 = 77.0 (^^7.0 = 13.7) 
and $73.4 = 167 (ilr.o = 92.3). 

FIG. 8. Width at half maximum Tg (diamonds) of the high-frequency resonance of the spectrum 
$''(aj) of the HSS for packing fraction ip = 0.60 as a function of the resonance-maximum position 
j^max £qj, various wave vectors. The straight line represents the small-wave-vector asymptote 
Tg = 7(0^"^/woo)^ = K"{uj = 0)(rj^"^)2 where K"{uj = 0) = 0.00182. 

FIG. 9. High-frequency sound-resonance position f2™"^ (diamonds) and resonance width F^ 
(circles) as a function of the packing fraction, determined for q = 1.8. The crosses exhibit the 
maximum position cop of the AOP of the density-fluctuation spectrum for wave vector q = 7.0. 

FIG. 10. Reactive parts of the moduli A{uj) = 1 — ujK'{uj) and spectra K"[lo) of the HSS for 
the packing fractions ip = 0.600 (solid), ip = 0.567 (n = 3, dashed) and ip = 0.540 (n = 4, dotted). 



FIG. 11. Mode-coupling coefficients Wg and Wg determining via Eqs. (10) the scattering 
and decay contributions, respectively, to the kernel K of the HSS glass states. The curves with 
label n = 3 and n = 4 refer to the packing fractions 99 = 0.567 and ip = 0.540, respectively. The 
insets show the coefficients for g < 3 magnified by a factor 100. 
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FIG. 12. The spectrum K"{lo) of the kernel for the HSS at packing fraction (p = 0.60 reproduced 
from Fig. |l^ (full line). The dashed and the dotted curves exhibit the one-mode contribution 
K^^'"{uj) and the two-mode contribution K^'^'"{u>), respectively. The full line with label R denotes 
the Rayleigh contribution Rquj"^ {Rq = 1.5 10^^^) magnified by a factor 10^. 



FIG. 13. The full line reproduces the spectrum K"{uj) from Fig. 12. The dashed line exhibits 
the stiff-glass-approximation result K"{lo) from Eqs. ( pO[) with Q = 120, wi = 0.388, u = 0. The 
dotted line is the extension of this approximation by incorporating v = VI? K^"^'" {u = 0) ~ 23. The 
dot-dashed line shows the extension with iv replaced by the kernel Vl? K^"^' {uj) from Eq. ([2l|). 

FIG. 14. Fluctuation spectra ^"[lo) for the one-component schematic model defined in Sec.IVA 
for a mode-coupling functional J-[f\ = vif + 1^2/^ (solid lines) and = 1. The states refer to 
the glass with distance parameters e = 1/4", n = 0,1,2 (compare text). The dashed line exhibit 
the stiff-glass approximation given by Eqs. ( pO[ ) with zv = 0. The dotted lines show the extended 
description including a z^ 7^ which was evaluated from the a; = limit of Eq. (pi]) . 



FIG. 15. Spectra ^^"{uj) of the solutions for the second correlator of the two-component model 
defined in Sec.IVA with t>* chosen so that f^ = 0.9. The first correlator needed as an input for 
the memory kernel ( pq) is the one corresponding to the n = curve in Fig. ^. The dashed lines 
are the stiff-glass approximations, Eqs. (^), to these curves. The arrows indicate, which axis 
corresponds to the curves. 

FIG. 16. Rescaled tagged-particle-density-fluctuation spectra 10^ ^g"{uj)/q'^ of the HSS for 
the packing fraction ip = 0.600 at wave numbers q = 3.4 (solid), q = 7.0 (dashed), and q = 10.6 
(dotted). 

FIG. 17. Density of states g{uj) of a HSS, calculated for the packing fraction 99 = 0.60. The 
dashed line shows 5(0;) = 9.2 x lO^^w^, describing the density-of-states asymptote at small fre- 
quencies. The arrow marks the Debye frequency ujd = 309 for the longitudinal sound. 
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FIG. 18. Dynamic structure factor S{q,u)) = Sq^g{io) of the HSS for packing fraction ip = 0.60 
as a function of the wave number q at some fixed frequencies uj. The Debye vector is qd = 4.08; 
threshold and maximum of the boson peak are located near uj = 45 and a; = 85 respectively 
(compare Fig. |l^). The lines are guides to the eye. 
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Figure 16: W. Gotze and M. R. Mayr 
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Figure 17: W. Gotze and M. R. Mayr 
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Figure 18: W. Gotze and M. R. Mayr 
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